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the hydrocarbon K 3 picene superconductor near the metal-insulator transition 



Minjae Kim and B. I. Min* 
Department of Physics, PCTP, Pohang University of Science and Technology, Pohang, 790-784, Korea 

Geunsik Lee, Hee Jae Kwon, Y. M. Rhee and Ji Hoon Shim* 
Department of Chemistry, Pohang University of Science and Technology, Pohang, 790-784, Korea 

(Dated: June 13, 2011) 

We have investigated the electronic structures and magnetic properties of a newly discovered 
hydrocarbon superconductor, Kapicene, having T c —18 K. We have shown that the metal-insulator 
transition (MIT) is driven in Kapicene by 5% volume enhancement with the formation of the local 
magnetic moment. Active bands for the superconductivity near the Fermi level (Ef) have hybridized 
character of LUMO and LUMO+1 of the picene molecule. Fermi surfaces of Kapicene manifest 
neither prominent nesting feature nor marked two-dimensional behavior. By estimating the ratio of 
the Coulomb interaction U and the bandwidth W of the active bands near Ef, we have demonstrated 
that Kapicene is located in the vicinity of the Mott transition. Our findings suggest that Kapicene 
is a strongly correlated electron system. 

PACS numbers: 74.20.Pq, 74.70.Kn, 74.70.Wz 



Recently, the first hydrocarbon superconductivity has 
been observed in the K-doped picene, K 3 picenei. Its 
transition temperature T c =18 K is comparable to that of 
the alkali-metal doped fullerene system (T c ^38K)2. In 
Kapicene, K atoms are intercalated in the stacked picene 
molecules, as shown in Fig. [Ha). A picene molecule 
(C22H14) consists of five connected benzene rings with 
armchair edge (see Fig.QJb)). Interestingly, alkali metal- 
doped Aapicenes (A=Na, K, Rb, and Cs) exhibit the 
cation A-dependent physical properties^. With increas- 
ing the cation size, the system varies from a Pauli-likc 
paramagnet for A=Na to superconductors for A—K and 
Rb, and then to an insulator for A=Cs. This feature sug- 
gests that Aspicene is susceptible to the chemical pres- 
sure. This trend is also reminiscent of A^Cqo (A=K, Rb, 
and Cs), in which changing cation from K to Cs results in 
6% volume expansion so as to induce the metal-insulator 
transition (MIT)£. 

Unconventional superconductors often manifest the in- 
cipient MIT and magnetic instability, which reflect the 
existence of strongly correlated electrons^. The par- 
ent compounds of cuprate high T c superconductors are 
typically antiferromagnetic (AFM) insulators^. Also the 
organic superconductors k-(BEDT-TTF) 2 X and A^Cqq 
(A—K and Cs) show the phase transition from an 
AFM insulator to a superconductor upon increasing 
pressure^. These examples suggest that investigation 
of the electronic structures of Kapicene under volume 
change would provide an important clue to the pairing 
mechanism of superconductivity. 

In this report, we have explored the volume-dependent 
electronic and magnetic properties of K 3 picene using the 
first principles band structure calculation within the den- 
sity functional theory. We have shown that the MIT 
and AFM transitions are driven in Kapicene by 5% vol- 
ume enhancement with the formation of the local mag- 
netic moment. Our results indicate that Kapicene itself 



is on the verge of both MIT and magnetic instability. 
We have found that the active bands for superconductiv- 
ity have hybridized character of LUMO and LUMO+1 
of the picene molecule. Quantitative estimation of U/W 
(£/:Coulomb interaction, W: band width) also provides 
that Kapicene is near to the Mott transition. 

In the unit cell of Kapicene, there are two stacked 
picene molecules and six K atoms (see Fig. (Ha)). The 
positions of K atoms are still undetermined experimen- 
tally. Hence, we first optimized the internal atomic posi- 
tions in Kapicene with retaining the observed monoclinic 
structure (lattice constants: a=5.912, b=8.707, c=12.97 
A, and a=92.77°)»i For this purpose, we have used the 
pseudo-potential band method of VASP package 8 . For 
the exchange-correlation potential, the generalized gra- 
dient approximation (GGA) of Perdew-Burke-Ernzerhof 
(PBE) was utilized^. We have used 100 k-points inside 
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FIG. 1: (Color online) (a) Crystal structure of Kapicene. 
There are two stacked picene molecules in a unit cell. K 
atoms are intercalated between the stacked picene molecules, 
(b) Molecular structure of a picene. 
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FIG. 2: (Color online) DOS for NM (blue solid lines) and 
AFM (red dotted lines) phases of Kapicene. Inset shows the 
amplified DOS near Ef- (a) Kapicene with the original vol- 
ume (KaP[OV]), and (b) Kapicene with 5% volume expansion 
(K 3 P[EV]). 
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FIG. 3: (Color online) (a) Band structure for NM phase of 
Kapicene with original volume (KaP[OV]). (b) Band struc- 
ture for AFM phase of Kapicene with 5% volume expansion 
(KaP[EV]). Notice that two bands near Ef are split into U 
and L which are represented by blue (dark) color. 



the first Brillouin zone. The convergence of the total 
energy with respect to number of k-points was checked 
to have precision of less than 4 meV per formula unit. 
For the structural optimization, the atomic relaxations 
are terminated when the forces in all atomic sites are 
less than 0.05 eV/A. The stable positions of K atoms 
are determined as four above the end benzene rings and 
two above the center rings with nearly two dimensional 
stacked arrangement of picene molecules^, similarly to 
the results of Ref. 11,12 . However, the arrangement of 
picene molecules in Ref/^i are perpendicular to ab plane 
losing the conventional herringbone structure of the pris- 
tine picene solid. Our optimized structure retains the 
conventional herringbone structure in accordance with 
RefJ^. We have also optimized the internal atomic po- 
sitions of Kapicene with 5% volume expansion in order 
to examine the chemical pressure effect, as is realized 
in Cs3picene. Then, based on the optimized structures, 
we have obtained the electronic structure of Kapicene by 
employing the more precise all electron full-potential lin- 
earized augmented plane wave (FLAPW) band method 13 
implemented in WIEN2k package^. To explore the mag- 
netic properties of Kapicene, we have considered the 
AFM spin configuration with the opposite spin polar- 
izations of two independent picene molecules in the unit 
cell. 

Figure [2] shows the density of states (DOS) of 
Kapicene. DOSs of K 3 piccnc with original volume (here- 
after K3P[OV]) and with +5% volume expansion (here- 
after K3P[EV]) are provided in their nonmagnetic (NM) 
and AFM phases. Pristine solid picene is a semicon- 
ductor with a wide band gap of 3.3 eV. 15 Figure [2(a) 
shows that K 3 P[OV] also has an energy gap structure 
(-1.7 eV) between HOMO and LUMO related bands. 



In Kapicene, three electrons are transferred from three 
intercalated K atoms into one picene, occupying LUMO 
and LUMO+I related bands, as was pointed out by Ko- 
sugi et. al. 11 . For K 3 P[EV] in Fig.rjjb), the band widths 
become narrower to have a larger energy gap (~ 1.9 eV) 
between HOMO and LUMO related bands. It is seen 
that, for both K 3 P[OV] and K 3 P[EV], the DOS at the 
Fermi level (Ep) is much reduced in the AFM phase with 
respect to that of the NM phase. Noteworthy is that 
K 3 P[EV] becomes an insulator in the AFM phase, hav- 
ing a small energy gap of ~ 0.03 eV. The NM and AFM 
phases are almost degenerate. Nevertheless, the mag- 
netic moment of one picene, ~0.35 /^b, is large enough 
to stabilize the magnetic ground state d 6 ' 17 These fea- 
tures suggest that K 3 picene (K 3 P[OV]) is on the verge of 
the MIT and magnetic instability. The AFM structure in 
K 3 P[EV] with the localized magnetic moment of a picene 
provides the explanation of the significant enhancement 
of M/H observed in Cs 3 picenei. Like other molecular 
solids (KO2, K-doped pentacene, Cs 3 Cgo) having local 
moment carrying 2p electron s 2 ' 18 ' 19 , the local moment in 
K 3 P[EV] arises from the unpaired electron in the local- 
ized molecular orbital. 

In Fig. [3fa) and (b), the NM band structure of 
K 3 P[OV] and the AFM band structure of K 3 P[EV] are 
plotted, respectively. The former exhibits a metallic state 
with two bands crossing Ep, while the latter shows an in- 
sulating state. The metallic band structure of K 3 P[OV] 
results from the almost degenerate U and L bands near 
Ep, while the insulating state in the AFM phase of 
K 3 P[EV] emerges due to the shift up and down of re- 
spective U and L bands that had become narrowed due 
to the volume expansion. Notice that the degenerate 
bands in the NM phase of K 3 P[OV] at k-points along 
the (100)-(101)-(110) directions are split in the AFM 
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FIG. 4: (Color online) Charge densities (CDs) of (a) U and 
(b) L bands in Fig. [3] CDs of (c) LUMO+1 and (d) LUMO 
orbitals of a picene molecule. CD of the U band shows mostly 
LUMO+1 character with small hybridization with LUMO, 
while CD of the L band shows hybridized character of LUMO 
and LUMO+1. 



phase of K3P [EV] . The overlap of tt-tt orbitals of adjacent 
molecules along the (100) direction produces the sym- 
metric hopping along this direction, resulting in the de- 
generacy on the (100) zone surface for the NM phases 20 . 
This degeneracy on the (100) zone surface is lifted by the 
AFM spin ordering for both K 3 P[OV] and K 3 P[EV], due 
to the broken symmetry along the (100) direction. This 
feature is reminiscent of the Slater-type MIT mechanism 
accompanied by the cell doubling, and might suggest that 
I<3P[EV] can be categorized into a band insulator. How- 
ever, similar feature was observed in the typical Mott 
insulating systems such as cuprate and nickel oxides. 
Therefore, in order to confirm that whether K3P[EV] is 
a band or Mott insulator, one needs to check the U/W 
value of a system, as will be discussed later. 

We have examined the orbital characters of U and L 
bands near Ep that are responsible for superconductivity. 
The charge densities (CDs) of U and L bands are plot- 
ted in Fig. Hta) and (b), respectively. For comparison, 
the CDs of LUMO+1 and LUMO of a picene molecule 
are also provided in Fig. HJc) and (d), respectively. The 
CDs in Fig. 0] reveal that U and L bands have the mixed 
character of both LUMO and LUMO+1, even though 
LUMO+1 character is dominant for U. The energy split- 
ting between LUMO and LUMO+1 of a picene molecule 
is known to be 0.15 eVi This gap size is comparable 
to the band widths of four conduction bands near Ep 
(U, L, and two lower bands below L in Fig. EJb)), which 
correspond to 0.21, 0.12, 0.12, and 0.13 eV, respectively. 
The comparable size of the gap and band widths indi- 
cates that the hybridization interaction between LUMO 
and LUMO+1 is large enough to yield the hybridized 
bands with mixed LUMO and LUMO+1 character near 
Ep. Thus, the active bands for superconductivity near 
E F have hybridized character of LUMO and LUMO+1 
molecular orbitals, in agreement with the earlier report 
by Kosugi et. aL— . This is contrary to the case of 
pentacene, in which the energy splitting between LUMO 



and LUMO+1 (~ 1.28 eV) is much larger than that of 
picene.— 

Fermi surfaces of K3picene are provided in Fig. El Note 
first that there is no prominent nesting feature in the 
Fermi surfaces of K 3 P[OV] (Fig. EJa)), which is in agree- 
ment with Ref Secondly, there is no Fermi surfaces of 
clear cylindrical shape along the c-axis, indicating that 
this system would not have marked two dimensional be- 
havior, despite its layered-like crystal structure. This 
feature is indeed consistent with the band structure in 
Fig. EJa), which shows apparent dispersions along the 
c-axis, (100)-(101) and (001)-(000). Figure E^b) and (c) 
illustrate explicit effects of AFM spin polarization and 
volume enhancement on the Fermi surfaces of K^picene, 
respectively. For the AFM phase of K 3 P[OV] (Fig. Efb)), 
the volume of Fermi surfaces is reduced much, with re- 
spect to that of the NM phase of K 3 P[OV] (Fig. EI>)). 
It is because of the reduced DOS near Ep accompanied 
by the AFM spin polarization. The volume enhancement 
also reduces the Fermi surface shown in Fig. Elc), 

due to the decreased band width. 

We have calculated electrical conductivity a using the 
Boltzmann transport theor y 21 ' 22 . The cr/r (r: scattering 
time) value of the NM phase of K3P[OV] is estimated to 
be about 0.75 (10 19 S/m- sec). The calculated anisotropy 
is only about two, which implies that this system is not 
so anisotropic, as addressed in Fig. El[a). The reduced 
Fermi surfaces due to the explicit effects of AFM spin 
polarization and volume enhancement result in signifi- 
cant drop in the conductivity. The cr/r values of the 
AFM phase of K 3 P[OV] and the NM phase of K 3 P[EV] 
turn out to be 0.25 and 0.13 (10 19 S/m- sec), respectively, 
which are only 33% and 17% of that of the NM phase of 
K 3 P[OV] (0.75xl0 19 S/m • sec). The conductivity anal- 
ysis indicates that both volume enhancement and AFM 
spin polarization are effective for the MIT. The volume 
enhancement leads to the heavy effective mass in conduc- 
tivity by reducing W, while the AFM spin polarization 
leads to the reduced DOS at Ep by splitting two bands 
around Ep farther. 

In general, the volume enhancement of solid decreases 
the band width W and increases the Coulomb repulsion 
U, resulting in larger U/W value. The increased U/W, 
when it becomes larger than the critical value (U/W) c , 
would drive the MIT through a Mott transition. In- 
deed, in Kipentacene, the formation of insulating state 
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FIG. 5: (Color online) Fermi surfaces of Kapicene. (a) NM 
phase of K 3 P[OV], (b) AFM phase of K 3 P[OV], and (c) NM 
phase of K 3 P[EV]. 
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is experimentally realized^ 3 -. We have estimated U /W of 
K 3 picene, using the method of Rcf. 2 - that was applied 
to the pentacene solid. The effective Coulomb interac- 
tion (U e ff) is given by subtracting the polarization en- 
ergy (Ep i) from the bare Coulomb interaction (Ubare), 
Ueff = Ubare — E po i. Ubare of a charged picene molecule 
is estimated from the second order M0ller-Plesset cal- 
culations with 6-31G* atomic orbital basis set for an 
isolated picene by considering different anionic electron 
occupations 2 ^^. E poi for a solid K 3 picene can be deter- 
mined by using the partial continuum model, in which a 
charged molecule is assumed to be surrounded by a cav- 
ity of homogeneous dielectric medium 2 ^. We have ob- 
tained the dielectric constant and the size of the cavity 
by fitting the calculated E po i to the experimental value 
of the pristine solid picene^l. In this way, we have got 
Ubare = 2.25 eV and E po i = 1.40 eV, respectively, and 
thus U e ff is estimated to be 0.85 eV. One can estimate 
the U/W for K 3 P[EV] by adopting W of L+U bands 
near E F (0.30 eV) for the NM phase of K 3 P[EV]. Then 
we have obtained U/W — 2.83, which is larger than 
(U/W) c = 1.732&3&, so that a Mott insulating state is 
realized in K 3 P[EV]. This large U/W value for K 3 P[EV] 
suggests that K 3 P[OV], which exhibits superconductiv- 
ity, would also be near the boundary of the Mott transi- 
tion. 

In a recent report on the K atom adsorbed picene 
molecule- 1 , it was argued that the charge transfer from 
K atoms to LUMO and LUMO+1 is not decisive due 
to the strong Coulomb correlation effect in picene. The 
above analysis for the Coulomb interaction, however, 
shows that U e ff in K 3 picene solid is only about 38% 
of Ubare in an isolated picene molecule. Therefore, dis- 
tinctly from the case of the picene molecule, the charge 
transfer from K atoms to LUMO and LUMO+1 related 
bands will be possible in solid picene due to the reduced 



effective Coulomb interaction. 

In conclusion, we have studied the electronic structures 
and magnetism of a first hydrocarbon superconductor, 
K 3 picene, by employing the first principles band struc- 
ture method within the density functional theory. Active 
bands for superconductivity are found to have LUMO 
and LUMO+1 hybridized character. Fermi surfaces of 
K 3 picene manifest neither prominent nesting feature nor 
marked two-dimensional behavior. The latter is also sup- 
ported by the rather small anisotropy in the estimated 
conductivity. We have demonstrated that K 3 picene is 
on the verge of the MIT and the magnetic instability. 
Therefore, one needs to consider a 3D model with strong 
Coulomb correlation effects to investigate the mechanism 
of superconductivity in K 3 picene. 

After submission of the manuscript, we became aware 
of a related work by G. Giovannetti and M. Capone, 
Ref.— , in which the AFM state is found to become sta- 
bilized in K 3 P[OV] by the inclusion of correlation effects. 
Also recent experiment confirms that the Coulomb in- 
teraction (U) in solid picene is about 0.85 eV, which is 
much larger than the band width, Ref. 32 . These results 
are consistent with the present results suggesting that 
K 3 picene is a strongly correlated electron system. 
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